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Abstract 

The thermal conductivity of a nanostructure is sensitive to its dimensions. A simple analytical 
scaling law that predicts how conductivity changes with the dimensions of the structure, how- 
ever, has not been developed. The lack of such a law is a hurdle in "phonon engineering" of 
many important applications. Here, we report an analytical scaling law for thermal conductivity 
of nanostructures as a function of their dimensions. We have verified the law using very large 
molecular dynamics simulations. 
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Thermal conductivity is a size-independent property for macroscopic scale materials, but 
becomes sensitive to sizes when the feature dimension is reduced to nano- or micro- meter 
scale 1 . This size-dependence is critical for nanostructure applications. For instance, it can 
be utilized to engineer high thermal conductivity and heat dissipation of microelectronic 
elements and thereby to effectively increase their density in a device^. It can also be utilized 
to engineer low thermal conductivity of thermoelectrics systems to improve energy conversion 
efficiency^. Previous workP^ has used a Matthiessen rule^ to relate thermal conductivity k 
to sample length L: 

11a 

k{L) ~ Kb + L ^ 
where is the thermal conductivity of bulk material and a is a size independent constant. 

While this rule has been successfully applied^, it is only applicable for heat conduction 
through the "thickness" L of a film (i.e., the sample is assumed to have an infinite cross 
section). It cannot be applied for heat conduction along a direction in the plane of the film, 
nor can it be applied for any nanostructures with more than one dimension at the nano- or 
micro- scale. Because a general equation for thermal conductivity of nanostructures has not 
been developed, some previous analysis has used the solution of Boltzmann partial differ- 
ential equations in order to explore the effect of nanostructure sizes^ffl^. This approach is 
complex and depends on an empirical estimate of the specularity of the free surfaces. It has 
not been applied to nanostructures with arbitrary dimensions in all three coordinate direc- 
tions, nor has it provided non-sectioned analytical solutions. The lack of a tractable scaling 
law has posed a hurdle in "phonon engineering" of many nanostructure applications. To 
overcome this problem, we have developed an analytical scaling law that explicitly expresses 
thermal conductivity of nanostructures as a function of dimensions in all three coordinate 
directions. We have also verified the law using large molecular dynamics (MD) simulations. 

Consider the unidirectional heat conduction through the length L of a box-shaped sample 
with a thickness t and a width W, as illustrated in Fig. [Ha). It is recognized that the 
size effect on thermal conductivity comes from the surface scattering of phonons, which 
diminishes as the distance from the surface is increased. Hence, we divide the sample into 
different regions with respect to the surfaces, Fig. [TJ First, the (y-z) cross section of the 
sample is divided into nine regions as shown in Fig. [ijb). These nine regions extend in the x- 
direction into nine small box-shaped pieces (referred to as pillars hereafter) as shown in Fig. 
pTa). The eight pillars surrounding the center pillar essentially form a shell whose thickness 
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is assumed to be d. Each pillar exhibits an apparent mean thermal conductivity throughout 
length L. At a fixed large value of d (say in the order of the phonon mean free path), the 
boundary environment of each pillar is independent of sample dimension t and W. As a 
result, the apparent thermal conductivity of each pillar is a function of L only. Fig. |T|(b) 
indicates that for isotropic materials, the nine pillars fall into three different types i — 0, 1, 2, 
where % refers to the number of the y- or z- surfaces bounding the pillar. Correspondingly, we 
have three distinguishable conductivity functions k (L), k\ (L), and k 2 (L). Because heat 
transports through the pillars in parallel, the overall thermal conductivity of the sample can 
be calculated as an area- weighted average leading to: 

k (t,W,L) 

= k (L) - [«o (L) - «! (L)] -(j + w) 

Ad 2 

+ [«o (L) + k 2 (L) - 2k, (L)} ■ — (2) 





FIG. 1: Heat conduction through the length L of a box-shaped material with a thickness t and a 
width W. 

Now we consider the thermal transport of each pillar. As shown in Fig. W[c), the pillar 
can be divided into a center section with a length of L — 2u and two end sections with a 
length of u). At a large given value of u, the two end surfaces do not interact. As a result, 
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the apparent thermal conductivities of the center and the end sections of the pillar are 
independent of the length L, and therefore can be represented respectively by two constants 
Kj C and Ki t£ , where subscript % is the pillar number, and c and e designate the center and end 
sections. Because heat transports through the center and end sections in serial, the overall 
thermal resistivity (inverse of thermal conductivity) of the pillar can be calculated as the 
length- weighted average resistivity: kJ 1 (L) = (1 — 2uj/L) ■ k^ c + (2u/L) ■ K^ e , which can be 
rewritten as 

* (L) = L r ^f (3) 



L + Si 

where Si = 2co ■ (/q iC — / n i>e is a positive constant. Substituting Eq. ^ into Eq. (|2]), 
we have an analytical scaling law: 

L ■ k c 



k (t,W,L) 

L ■ K,Q,c 



L + 5 

L ■ Ki, c 



+ 



L + 5 L + 5 l 

L ■ Kq,c L ■ K2 yC 



2 • 



2d 2d 

T + w 



Ad 2 
t-W 



r r - , r,- ( 4 ) 

L + 5q L + 5 2 L + 5 l 
Eq. Q is valid if a sufficiently large d is used to subsume the surface scattering affected 
region. Once d is given, Eq. Q involves six parameters k 0iC , k 1c , « 2 ,o <$u $2i where k c 
is essentially the bulk thermal conductivity Kb, «i, c is the thermal conductivity near a flat 
surface, and k 2 , c is the conductivity near a corner region, see Figs. [T[b) and|l](c). For MD 
applications, it is necessary to perform several simulations at different dimensions in order 
to fit Eq. Q. If the minimum dimensions used in these simulations are represented by 
t m i n and W m i n , then the largest d that still enables all the MD data to satisfy the geometry 
condition (i.e., < 2d < t and < 2d < W) is t min /2 or W min /2 depending on which 
dimension is smaller. 

The model concepts described above can be applied to any sample shapes. Eq. (pj also 
has more general uses. For instance, we found that substituting t = W = 2r, K\ )C = «2, c , 
and 5i = 5 2 into Eq. Q resulted in the same axial thermal conductivity of a circular wire as 
a function of wire radius r and length L as we would otherwise derive by directly applying 
the concepts to the wire case. 

When t — 7- oo and W — > oo, Eq. (|4]) indicates that the inverse of thermal conductivity 
along the length of sample with an infinite cross section is a linear function of the inverse of 
the length 1/L, exactly matching the established equation, Eq. (TyHHS. When L — > oo and 



t — > oo or L — > oo and W — > oo, Fig. [T] corresponds to heat transport along a direction in 
the plane of a large film. Eq. Q then shows that thermal conductivity is a linear function 
of the inverse of the film thickness (W and t in the two cases both correspond to film 
thickness). Eq. Q can be effectively verified by checking these linear relationships using 
either experiments or MD simulations. Here we perform "direct method" MD simulations^ 
to verify the linear relationships. The [0001] thermal conductivity of a wurtzite GaN crystal 
was calculated at a temperature of 300 K. GaN is of interest because many of its applications, 
such as laser diodes and high electron mobility transistors^^, operate at high current and 
power densities where heat dissipation is crucial. 

A Stillinger- Weber potential developed by Bere and SerraP^^ was used. The computa- 
tional cell is aligned so that the x-, y-, and z- coordinates correspond, respectively, to [0001], 
[1100], and [1120] directions. A periodic boundary condition was used along the z- axis to 
simulate an infinite width of W — > oo, and a free boundary conditions is used in the y- direc- 
tion to simulate the commonly encountered [1100] surface. Series of thermal conductivities 
at two lengths (in the x- direction) of L = 260 A and 390 A and different thicknesses (in the 
y- direction) of t between 276 and 829 A and t — > oo (periodic boundary conditions) were 
calculated. All the simulations applied a very long averaging time of at least 11 ns (some 
reached 21 ns). Our systems are also relatively large (up to 900,000 atoms). Both a long 
averaging time and a large number of atoms available for averaging help generate highly 
accurate results^ that strengthened the conclusions. The calculated values of k are shown 
in Fig. [2^a) against 1/t, where the lines are produced using Eq. Q with the assigned value 
d = 138.13 A and the fitted parameters k , c = 178.38 W/K ■ m, k x . c = 151.65 W/K ■ m, 
5q = 1288.10 A, and 5\ = 1329.65 A. It can be seen that the MD data well satisfied the 
predicted linear relationship and the agreement between the MD data and Eq. Q is excel- 
lent. Previously calculated thermal conductivities at different sample lengths L but a fixed 
sample width W — >■ oo and a fixed sample thickness t — > oc^ are reproduced in Fig. [2^b) 
using the 1/k vs. 1/L scale, along with the line generated using Eq. Q with the same 
parameters. Again the predicted linear relationship is well satisfied and excellent agreement 
is obtained with only one set of parameters (d, kq^ c , ki jC , 5q, and 5i) for both thickness and 
length functions. 

In summary, we have developed an analytical law for size effects on thermal conductivity 
of nanostructures. This law is very well verified by MD simulations. We expect that it will 
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FIG. 2: GaN thermal conductivity as a function of sample dimension. 

enable fundamental methods such as MD simulations to be used to study thermal transport 
at realistic length scales, which would be otherwise impossible due to the limitation of the 
length scales that can be directly simulated. We also expect that this law can guide ex- 
periments to design nanostructured thermal devices. Note that when experimental thermal 
conductivity data is obtained at different dimensions, the same approach can be used to fit 
Eq. Q. Since experiments are likely to be performed at larger sample dimensions, larger 
values of d can be chosen to produce even more accurate results. Nonetheless, simulations 
performed here strongly indicated that even d = 138.13 A is sufficient for GaN. 
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